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We derive a common mathematical formulation for the eigenfunction statistics of Hermitian op- 
erators, represented by a multi-parametric probability density. The system-information in the for- 
mulation enters through two parameters only, namely, system size and the complexity parameter, a 
function of all system parameter including size. The behavior is contrary to the eigenvalue statistics 
which is sensitive to complexity parameter only and shows a single parametric scaling. 

The existence of a mathematical formulation, of both eigenfunctions and eigenvalues, common to 
a wide range of complex systems indicates the possibility of a similar formulation for many physical 
properties. This also suggests the possibility to classify them in various universality classes defined 
by complexity parameter. 
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I. INTRODUCTION 

The eigenfunction correlations of various generators of dynamics contain a wealth of information about the system 
e.g. localized or delocalized nature of the dynamics, decay rate etc. Recently the correlations were shown to be 
relevant for description of fluctuations of physical properties e.g. conductance in mesoscopic systems, peak-height 
statistics in coulomb blockade regime of quantum dots [l], Q • The correlations may vary from level to level or fluctuate 
in different realizations of a complex system. The strong fluctuations of eigenfunctions are already known to be the 
hallmark of many critical point studies e.g. metal-insulator transition in disorder systems [|[, spin glass and stock 
market fluctuations Q etc. Recent studies have revealed existence of the fluctuations in a wider range of complex 
systems e.g. in the area of quantum information, nanotechnologyQ and complex networks etc. As a consequence, 
a detailed information about the eigenfunction statistics of complex systems is very important and desirable. 

During recent years, many attempts have been made to statistically formulate the eigenfunction correlations of 
complex systems; see, for example, [H, i, 0, 0, ©, Gil El El El El El El El [3, [II]. 0ne of the main tools used 
in this context is the random matrix approach which can briefly be described as follows (see [§} for details). The 
presence of complicated interactions (among its various sub-units) in the system under investigation often makes it 
impossible to exactly determine the relevant operator e.g. in a matrix representation. The elements (some or all) 
of the operator in the representation can then be best described by a probability distribution. This permits one to 
replace the operator by an ensemble of the operators which is supposed to describe the generic properties and is 
referred as random matrix model of the operator. (In this paper, we focus on systems where complicated interactions, 
of any origin, lead to a partial or full randomization, of the operator, thus, allowing one to use a random matrix 
approach) . 

The choice of an appropriate random matrix model for a system is sensitive to its physical conditions (i.e nature 
and degree of interactions in various sub-units, symmetry and topological conditions, dimensionality etc.). This 
is because the distribution parameters of each matrix element depend on the interactions between related basis 
states (or parts of the system) which in turn are governed by system-conditions. In past, this has motivated an 
introduction of a variety of random matrix ensembles as models for a wide range of complex systems e.g. nuclei, 
atoms, molecules, disordered and chaotic systems, quantum chromodynamics, elastomechanics, electrodynamics (see 
reviews 0, i !, HI H HI, HI H, E3 and references therein for details), mathematical areas such as Riemann 
zeta function, enumeration problems in geometry and fluctuations in random permutations [251 ]. biological systems 
[HI, stock markets @, atmospheric sciences etc [3(j. For example, the systems with delocalized wave dynamics 
(extended throughout the system) and antiunitary symmetries can be well modeled by Wigner-Dyson ensembles; the 
latter are the Hermitian ensembles with Gaussian distributed matrix elements, with ratio of diagonal to off-diagonal 
variance a = 2, (originally introduced by Cartan [3lj . later developed by Wigner and Dyson to model compound 
nuclei and other systems) 8, 9]. The cases with partially violated anti- unitary symmetries can be well-modeled by 
Dyson's Brownian ensembles (BE) [1, HI] (see section VI also). The ensembles with arbitrary a 2) [33[, banded 
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matrices |23l . |34| (elements with non-zero variance within a band around main diagonal) and sparse matrices (with 
many elements with zero variance) have been successfully used to model statistical properties of the energy levels 
and eigenfunctions of systems with localized wave dynamics (e.g. quasi one dimensional wires and disordered systems 
of higher dimensions, chiral systems) [H, Hoi . [20L l35l . [36J. During last decade, many new ensembles were introduced 
to model the systems with unitary symmetries e.g. block form matrices for the cases with parity violation and pre- 
compound nuclei, chiral ensembles for systems with chiral symmetry in quantum chromodynamics 0, [26|], C and CI 
ensembles for cases withparticle/hole structures [£ j_ 37 [ , superconductivity etc [37j] . The non-Hermitian operators 
e.g. scattering matrices PI J$L [27| . transfer matrices 27j| or correlation matrices (appearing in time-series analysis e.g. 
stock market brain (28|, atmospheric studies [2£|) can similarly be modeled by circular ensembles [H, Ginibre 
ensembles (38[ and their more generic forms [39l |4C |. (The breadth of the subject is such that it is not possible to 
give a detailed account of all ensembles or include all references here). 

The applicability of random matrix ensembles to complex systems has been under investigation for past few decades. 
The validity of the models, however, has been extensively verified in context of the eigenvalue fluctuation only; see 
reviews 0, Q, !, [H, H& H [H, EH, S3, S3 and references therein. The validity in the domain of eigenfunction 
fluctuations is so far mostly studied either in ergodic regime of the wavefunctions (see (ll. li~2l. [i~3l. [IH liH li~7[ [l8LJl9l. HH 
for some original papers and reviews) or for quasi 1-d systems [l(| and specific cases e.g. disordered systems^fTThe 
growing technological demands as well as the observations of hitherto unknown features among eigenfunctions (e.g . 
multifractal structures at critical points) of a wide range of complex systems (for example, see [H, 0, HH, |45| . |46| |) 
have made it imperative to seek the statistical information in higher dimensions and beyond ergodic regime. This 
motivates us to pursue the present study. It is also desirable to explore the possibilities of any connection among 
the critical point behavior of the eigenfunctions of different complex systems. One way to show the connection i s by 
describing their various measures by a common mathematical formulation if possible. A recent study [49|, IHol l5l| . 
has shown the existence of a similar formulation for the case of level-statistics where system information enters 
through a single parameter, basically a function of all system parameters. The well-known connection between the 
statistics of eigenfunctions and eigenvalues in non-ergodic regime [3[ motivates us to seek a similar formulation for 
the eigenfunctions too. Such a formulation can also be useful in deriving the measure of one complex system from 
another. 

The paper is organized as follows. The section II contains a brief revision of the single parametric formulation of 
the multi-parametric probability density of matrix elements for a wide range of complex systems (see (49j for details) . 
The section III describes the derivation of the complexity parameter governed diffusion equation for the eigenfunction 
components (of the same eigenfunction or different ones) which is used in section IV to study the distribution of some 
of the important fluctuation measures. The other measures can also be derived following the same route. Although, 
the diffusion approach seems to complicate the calculation by introducing a dependence on the initial conditions, 
however, as discussed in section IV, the statistics of the system can be recovered by integrating over all physically 
allowed initial conditions. The approach has an extra advantage: it provides a common analytical base for the 
systems which can be modeled by our ensemble (given below by eq.(l,2)) The section IV briefly discusses the role of 
complexity parameter in various transitions induced due to change in system-specifics. The section V contains details 
of the numerical verification of our analytical claims. We conclude in section VI by summarizing our main results and 
their potential applications. 



II. SINGLE PARAMETRIC FORMULATION OF THE MATRIX ELEMENTS PROBABILITY DENSITY 

The eigenvalues and eigenfunctions of an operator, say H, of a system can be obtained by solving the eigenvalue 
equation HUi = X JJi (with Ui and \ as the eigenfunction and corresponding eigenvalue respectively) and any other 
physical information can then be deduced, in principle, from this knowledge. In the case of a complex system, however, 
the exact form of an operator e.g Hamiltonian may not be known or it may be far too complicated to solve. To deal 
with such a situation, one has to make a statistical hypotheses, known as maximum entropy hypothesis, for H [47| : 
a sufficiently complicated system can be described by a matrix which is as random as possible under the conditions 
compatible with the nature of the dynamics as well as the symmetry requirements. Thus if the symmetries and the 
nature of the operation is approximately known in a basis space preserving the symmetries, it can be modeled by an 
ensemble of full or sparse random matrices in that basis. For example, an equal probability of dynamics in each region 
of a specific space suggests a uniform spread of the eigenfunctions in the entire associated basis space. The operator 
in such a basis will therefore be a full matrix, (k\H\l) = XiUkiUu being of the same order for all combinations of 
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basis vectors \k), \l) (with Uf-i as the k th component of eigenvector [/$). On the other hand, the dynamics localized 
in a space leads to variation of the eigenfunction intensities in the associated basis and the operator will be a sparse 
matrix. 

It is clear from the above that, contrary to eigenvalues, the eigenfunction statistics depends on the basis in which 
the matrix is represented. The knowledge however is still relevant because (i) it can provide important information 
about the system dynamics in a given basis-space of interest, (ii) it is also possible to define a relevant basis to 
represent an operator: it is the basis in which the constraints on the operator appear in a natural way. For example, 
for the time-reversal invariant systems with integer angular momentum, the relevant basis is the one in which their 
Hamiltonians are simultaneously expressed as real-symmetric matrices |24| . 

In this paper, we consider a prototype distribution which can model a wide range of complex systems, namely, an 
ensemble of Hermitian matrices H, each of size N, described by a Gaussian probability density 



N 

P (H,h,b) = Cexp[-Y / J2 (l/2h H ., s ).(H kl . s - b kl;s ) 2 } (1) 

s=l k,l=l;k<l 

Here \k), \l) are unit vectors of the arbitrary basis of size N, chosen to represent H with H k i = (k\H\l). The subscript 
"s" refers to the components of Hki, (3 as their total number ((3 = 1 for real variable, (3 = 2 for the complex one), 
C as the normalization constant, h as the variance matrix with hki }S — (H k i s ) and b as the mean value matrix 
with (Hki-, s ) = bki-s- Our choice of Hermitian nature of the ensemble restricts the present discussion to the class of 
systems with conservative dynamics. Following maximum entropy hypothesis, the above ensemble can well-describe 
the distribution of the operators for which the average behavior of the matrix elements and their variances is known. 
Based on the complexity of the system, the elements of the parametric matrices h, b can have various functional forms 
e.g. exponential, power law etc. For example, the limit hkiys ~~ > 0, corresponds to non-random nature of Hki- The 
limit hki; S — * 7 1 bki is — > for all {k, l,s} gives the density for a Wigner-Dyson ensemble [![: p{H) oc e~ Trff2 . The 
limit hki- t s —> [a-dfiki + ct (l — $ki)}, bki- s —> for all {k,l,s} gives the density for a Rosenzweig-Porter ensemble [331; 
(note Brownian ensembles have the same density-form too, see section VI). A band matrix ensemble (lOL [20L l23l l34j 
with Gaussian distributed matrix elements and band length t can be obtained by substituting hki- a — ► 0, bki- s — ► if 
\k — l\ > t and hf-i-s oc a(\k — l\) for \k — l\ < t with various possible forms of function a e.g. exponential, rectangular 
etc. Similarly other ensembles with uncorrelated matrix elements, some of them with Gaussian randomness and others 
non-random, can be represented by appropriate choice of h and b parameters [50j . 

The eq.(l) is applicable for the cases with mutually independent matrix elements, with no condition imposed on 
their moments higher than 2 nd order. Here we briefly mention only two such cases, namely, disordered systems and 
mixed dynamical systems; (the application to other cases e.g. algebraic or algorithmic complexity [55| will be discussed 
elsewhere). During recent past, specific cases of eq.(l) have been extensively used to model the energy level statistics 
of disordered system within independent electron approximation; (the latter results in independent matrix elements 
of the Hamiltonian) . One such example is the power law random banded matrix (PRBM) ensemble (each hki with 
a power law dependence on the distance from the diagonal) [3~i| which has been shown to be a good model for the 
level-statistics of Anderson Hamiltonian (AH) 3]. The ensemble (1) was also used recently to prove, analytically as 
well as numerically, the single parameter scaling of the level-statistics of Anderson Hamiltonian and its mapping to 
single parametric Brownian ensembles [5(| ■ (Later, in this paper, AH, BE and PRBM cases are also used to verify 
our analytical predictions for the eigenfunction statistics). 

Another potential application of eq.(l) is to systems with mixed dynamics where, similar to disorder, KAM tori lead 
to a localization of dynamics (48|. The connection of quantum systems in classically chaotic and inte grab le regime 
to Wigner-Dyson ensembles and Poisson ensembles, respectively, is already well-established 0, El, 52, 53]. In past, 
it has been suggested that a mixed Hamiltonian (or time-evolution operator) in a relevant basis should appear as a 
block diagonal matrix, each block being associated with an isolated region of the classical phase space [2J, [52|, H3] ■ 
In cases, where a chaotic region can be decomposed in nearly, but not completely isolated subregions, blocks are 
expected to be connected through small but non-zero matrix elements. The average size of these matrix elements i.e 
the quantum constraints will be related to flux connecting different regions i.e to classical information. We further 
suggest that the regimes with stable islands can be modeled by the blocks with non-random elements (e.g. zero 
variance and non-zero mean). The chaotic regimes can be modeled by blocks with randomly distributed elements 
(e.g. same non-zero variance for all elements or only within a band). Based on nature of the dynamics, a chaotic block 
may further contain hierarchy of random and non-random sub-blocks; various diagonal blocks may also be correlated. 
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Eq.(l) can then be applied to model the Hamiltonian by choosing the matrix element variances appropriate to the 
block in which they appear. 

Eq.(l) can not serve as a good model for the cases with correlated matrix elements. For example, particle-particle 
interactions in nuclei [l|, 0, HJ, H3| and electron-electron interaction in disordered systems can lead to correlations 
among elements of the Hamiltonian [l|, l5ll ]: the correlation coefficients depend on various system parameters. In 
general, such cases can occur when the interaction (described by H) between any two basis states is influenced by 
the other states. In past, consideration of particle correlations in nuclei led to introduction of embedded ensembles 
@)SHIi; however no significant progress has been made so far in dealing analytically with these ensembles. 

In general, an increase of constraints on the system-dynamics subjects higher moments of the matrix elements to 
certain specific conditions. This motivates us to consider an alternative ensemble, namely, maximum entropy ensemble 
with restricted higher moments. Within maximum entropy hypothesis, the probability density for such cases turns 
out to be non-Gaussian [5l|: p(H) = Cp{H) where 



n 
s—l r—1 



-J2 b P( r ) (II, , h 

p(r) 



(2) 



with C as a normalization constant. Here each Hjk is expressed in terms of its (3 components, (f3 = 1 for the real- 
symmetric matrices and (3 = 2 for the complex Hermitian case): Hjk = 'J2s=i(i) s ~ 1 Hjk; S - Here symbol p(r) refers to 
a combination of r matrix elements chosen from a set of total M = N(N + l)/2 of them; note the terms present in a 
given combination need not be all different. The TT r ■ . implies a product over r terms present in the p th combination 

with coefficient b p (r) as a measure of their correlation: i i ^i P o P ;a) = gjrTT - ^he J2 v (r) * s a sum over au possible 



combinations (total (M) r ) of r elements chosen from a set of total M = N(N + l)/2 of them. 

The potential use of eq.(2) to disordered systems with e-e interaction is discussed in [5l| . Here we briefly discuss 
few more examples. The systems with chiral symmetry can be modeled by Hermitian ensembles with block form 

{^W + ^0^)' W ^ ^ as a m& trix of size N. Here, due to H^\i = H* k , for 1 > k,l > N (with 
I + N), the correlations between these elements are subjected to condition: 



matrices H = 
k' = k + NA'i 



- (\H k ',i\ 
Hk,i) = (Hk',i>) 



\H^\ 2 ) 



(3) 



However, due to H^.i = Hk'.i' = 0, all other matrix elements are uncorrelated. For a simple explanation, let us restrict 
to a case of real matrix W with Gaussian distributed elements. The ensemble can then be represented by eq.(2) with 
ri = 2 or cquivalently by density, 



p{H, a, b) = Cexp 



"■■< u 

i<j,k<l 



HijHki — ^ akiHk 



ki 



(4) 



with following conditions on a and b: 



bk',t,i,k' = bk',i,k',i = bi,k>,i,k>, 

au = a k',i' = 0, o>k',l = ai k' 



(5) 



Note, b parameters corresponding to other pairs of elements (both or at least one in diagonal blocks) diverge due 
to zero correlation between elements in such pairs. The cases with other types of correlated blocks can similarly be 
modeled by applying appropriate conditions on the b parameters which correspond to combinations of matrix elements 
appearing in opposite blocks. For example, the ensemble C describes the cases with particle-hole symmetry with a 

i / | ([S3] for details). Now the correlations between various elements must be subjected 



Hamiltonian H 



B 



to conditions Hk',i = H, 
eq.(3)): (H kll H v #) 



-A 

l,k>> H Kl 
<\HkA 2 ) 



-Hii y k r , this implies another set of non-zero correlations (besides those given by 
-(\Hi' y k'\ 2 )- For Gaussian distributed real matrices A and B, the case can again 
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be modeled by eq.(4) however now b parameters for other pairs (besides those given in eq.(5)) can also be finite and 



J k.l,k,l 



■>2',fc',2',fc'i 



(ii)b k 



Lk',1 



?k,l,l,k' 



alphab, 



V,k',l,k' 



ab 



V,k>,k',l 



with a 



1. 



satisfy the equality: (i)abk.i.i 
Recently it was shown [4£ 

stages of a Brownian type diffusion process in the matrix-space, evolving with respect to a single parameter which is 
a function of the distribution parameters of the ensemble: ... 



51j that the distribution p for both cases (eqs.(l,2)) appear as the non-equilibrium 



dp_ 

BY 



with 



^ ± = E 



k.l:. 



= L + p 

9ki 9 
2 dH kl , s 



± lH kl , s 



(6) 



(7) 



where g k i = 1 + Ski- The variable Y is the parameter governing the evolution of matrix elements subjected to various 
system conditions. For case (1), 



Y 



2M 7 



-In 



k<l s=l 



\bki-s? 



C 



(8) 



with Y[ implying a product over non-zero bki-, s and xu-s = 1 — (2 — 5ki)jhki;s, Co as a constant determined by the 
initial distribution, M as the number of all non-zero parameters x k i- s and b k i- s - The parameter 7 is arbitrary, giving 
the freedom to choose the end of the evolution; the Imihki-s — > 7, b k i- s — * for all fc,^ gives 1" — > 00 and the steady 
state (a Wigner-Dyson ensemble). The distribution parameters being indicators of the complexity of the system, Y 
can be termed as the complexity parameter [49| . Some examples of the calculation of Y from eq. (8) are discussed in 
[49l | (for banded ensembles) and in (50| (for Anderson Hamiltonian) . The Y in case of a mixed system can similarly 
be calculated if one knows the details of the mixed dynamics. 

In general, the form of parameter Y for eq.(2) is quite complicated; its details can be found in [5l|. However, for 
case (4), that is the Gaussian version of case (2), Y can be given as 



Y 



E 

kl 



E 

ijkl 



db 



ijkl 



X + constant 



(9) 



where summation is implied over the distribution parameters with finite values only, and, X = \Y] k i /fc; +X)yfe; fijkl) 1 
with fu = ja k i - 2[Y, mn a mn b k i mn + a mn bkinm] , fijkl = hhjkl ~ 2 J2 mn kjnmhlmn] ■ For further clarification we refer 
the reader to [51( where an example, namely, the modelling of lowest Landau level of a disordered quantum Hall 
system by eq.(4) and calculation of corresponding Y is discussed. 

It is easy to solve eq.(6) for arbitrary initial condition, say Hq at Y = Yq,: p(H,Y\Hq,Yq) cx exp[— (a/2)Tr(H — 
rjHo) 2 ] with a = 7(1 — ?/ 2 ) -1 and rj — e lY . The probability density of H can now be extracted by integrating over an 
ensemble of initial conditions: p(H, Y — Yq) = J p(H, Y\Ho, Yo)p(Hq, YojdHo. It is often useful to study the statistics 
of the perturbed Hamiltonian H in the eigcnfunction basis of unperturbed Hamiltonian Hq. Thus if the eigenfunctions 
of Hq are chosen as the basis vectors \l) etc, and, the initial distribution is given by p(Hq) oc exp[— (1/2) Y^j 
the eigenvalue equation UH — AU can be used to transform p(H) from matrix space to eigenvalue-eigenvector space 

{A.r|, 



P(H,Y) 



N 

n 

k,l;k<l 



Xifexp 



N 

-(7/2) 



(7M/2)El A fc- A '| 2 |^| 2 | [/ : 



k<l 



(10) 



where p = [e 2 ^-^ - l]" 1 . 

As indicated by eqs.(6) and (10), the ensemble densities for various complex systems (i.e. different h,b matrices) 
undergo a similar evolution as a function of Y. The Y"-governed flow for the joint distribution of the desired eigen- 
functions components and eigenvalues can be obtained, in principle, by integrating either eq.(6) or eq.(10) over all the 
undesired ones; however it is easier to integrate eq.(6). To explain the technique, we consider some of the important 
cases in this paper. 
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III. DIFFUSION EQUATION FOR EIGENFUNCTION COMPONENTS AND RELATED 

EIGENVALUES 



The k component Uu of an eigenstate Ui is a measure of the contribution of k basis state to the eigenstatc. 
Experimental observations of complex systems indicate the level to level variations as well as sample to sample 
fluctuations of the contribution. As a result, the knowledge of average behavior of the components is not enough and 
one needs to study their distribution. In this section, we consider the joint probability distributions of a few relevant 
combinations of the components of the operator H . The basis chosen for the representation of the eigenfunctions 
is the one in which the matrix elements of H have distribution (1) (or (2)). We use following notation in reference 
to various correlations: For a joint distribution P rs , the subscripts r and s refer to the number of components of 
each eigenvector and the number of eigenvectors considered, respectively. For example, for a joint distribution of n 
components of m eigenvectors along with their eigenvalues, r = n and s = m. 



A. Joint Distribution of a Given Component of All Eigenfunctions and Eigenvalues 



It is often relevant to know the influence of a particular basis state on the system dynamics at various energies 
and with varying complexity of the system. The information can be obtained by a knowledge of the distribution 
of the same component of various eigenfunctions and its Y governed evolution. For example, let us calculate the 
joint distribution of a given component of all eigenvectors and the eigenvalues. Let Pin(Z, E, Y) be the probability, 
at a given Y, of finding the j component Uj n of the eigenfunctions U n of H between Zj n and Zj n + dzj n and the 
eigenvalues A„ between e„ and e„ + de„ for n = 1 — > N (with Z = {zj n },E = {e„}) . It can be expressed as an 
average over entire ensemble p: 



P 1N (Z, E,Y)= J f N (Z, E, U, A) p(H, Y) AH (11) 



with f r (Z, E, U, A) = n«=i $( z jn — U 3 n)S l3 ~ 1 (z* n — U* n )S(e n — A„). The Y dependent evolution equation for P iN 
can now be derived by connecting the parametric derivatives of Pin to its derivatives with respect to eigenvectors. 
The steps can briefly be described as follows: As F-dependence of Pi at comes only through p, one can write 

= J dHf N L+p = J dHp L-f N + jPin (12) 

with 7 = f3N(N + 1)7/2. Eq.(12) is obtained, first, by differentiating eq.(ll) with respect to Y, then using eq.(6), 
followed by partial integration. Due to 6- function nature of Jn, its derivatives with respect to matrix elements can 
further be reduced to the derivatives with respect to Z and E, 



Of. 



N 



dH k i: S 



N 

E 



d\ n df N t dU jn df N t dU* n <9/_ 



N 



dHki-s de n dHki-s dz-. 



(13) 



The 2 nd derivative of f N can now be obtained from eq.(13) (see [49j]). The substitution of eq.(13) in eq.(12) helps as 
the derivatives with respect to Zj n and e ra can be taken out of the integral. It can further be simplified by a knowledge 
of the effect of a small perturbation of H on its eigenvalues and eigenvectors; the related results are given in Appendix 
B (see [49[ for the details). Using the relations, eq.(12) can be rewritten as 



dP 

— = (L Z + L* Z )P + L E P (14) 



where P = C\P\n, C\ = e lY , and 
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L, = 



ft 2 



N 

E 



n,m—l;n^m J 



<9e„ 



7e« + 



. i-jm 
J" 



<9e„ 



(15) 



where implies the complex conjugate of L^; note = f° r = 1 case. Eq.(14) describes the F-governed 
diffusion of a given component of all eigenvectors and all eigenvalues. Its solution depends on the choice of initial 

— (1/2) \ ' H 2 

condition H a . In the diagonal representation of H (taking p(Hq) oc e 0;iJ ', the solution can be given as 



N 



Pin oc J| \e k - e^exp 



k,l;k<l 



N 



-(1/2) £e 2 - W2) E l e « - e m | 2 |^ n | 2 |^-. 



(16) 



with /z same as in eq.(10). (Note, the above result can directly be obtained from eq.(10) too). 



B. Joint Distribution of all Components of A given Eigenfunction and its Eigenvalue 



The distribution of the components of a specific eigenstate contains information about various basis states con- 
tributing to the state which in turn determines its spread. Proceeding along the same lines as for Pin, the diffusion 
equation for the joint probability Pjvi of the components U nk , n = 1 — > N, of an eigenvector U k and the corresponding 
eigenvalue can also be obtained. The evolution of 



P N1 (Z k ,e k ,Y) = J f kP (H,Y) dH, 
with f k — 5(Zk — Uk)5^^ 1 {Zl — U k )S(e k — X k ), can again be shown to be described by 

dP N i 



dY 



F k + F* k + L ek P N1 



where F k = (f3 2 /4) £ 2 =1 L qk with 



N 



L\ k 

L2k 
Le k PNl 



d 



E a-y [ Znk Qnn;k\ 
"Z nk 



and 



Q 



mn\k 



n=l 
N 

E 

m.n—l 

d 

de k 



E 



d' 2 



dz nk dz* mk 

dP 

de k 



Q\1 



l^P + + 0Q°nn;k 



( z njZ m j) T 



PN2 dTj. 



(17) 



(18) 



(19) 



(20) 



Here d Tj = dejd^Zj with d Si Z = dZdZ* and Pn2 = PN2(Z k , Zj,e k ,ej) is the joint probability of all components of 
the two eigenvectors Zj = {z n j} and Z k = {z nk } (n = 1 — > N) alongwith their eigenvalues ej and e k , respectively: 



Pn2 = J f k f jP (H,Y) dH. 



(21) 
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The presence of eigenvalue-eigenfunction correlations in exponent of p(H) (e.g. eq.(10)) as well as the terms of type 
(ej — e k )~ 2 in the denominator of eq.(20) makes it difficult to write F k (in eq.(18)) as a function of P^i{Z k , e k ). To 
write eq.(18) in a closed form, it is necessary to approximate Q™ n . k (Appendix A): 



Q 



mn:k 



D^(N-iy- r x s/2 [6r, 



'"ink 



^nkYP 



Nl 



(22) 



with D k as the local mean level spacing at energy e k . Here x = 1 f° r A* < C k > X ~ (/VOt) f° r A* > Cfc with £ as the 
ensemble averaged localization length of the eigenfunction U k and d as the system-dimension. The length £ enters in 
the formulation due to its relation with typical intensity of a wavefunction: |C/„fc| 2 = \z n k\ 2 ~ Cfc 

The substitution of eq.(22) in eq.(19) helps to write F k in terms of Pni, thus reducing the evolution equation (18) 
for Pjvi in a closed form: 



F k = 



AD 2 



N 

E 



d 



dz nk 



E 



dh 2 

dz, 



mk 



(23) 



with hi = (N - l)xz n kPNi, h 2 = x{&mn - z nk z* mk )P N i 



C. Joint Distribution of all components of q Eigenfunctions and their eigenvalues 

For certain physical properties e.g susceptibility, a knowledge of the correlations among two (or more) eigenvectors 
at two different space points may be required. The fluctuations of such correlations can be determined by the joint 
probability density P^ q of the components U nk (n = 1 — > N) of q eigenvectors U k (k = 1 — ► q) where 

r q 

P Nq (Z u Z 2 ,..Z q ,Y)= ]]f k p(H,Y) dH, (24) 
J k=i 

Proceeding exactly as in previous two cases, the y-governed diffusion of P/v g can be shown to be described as 



dP, 



Nq 



dY 



= Y,[h+F* k +L 



e k PNq 



k=l 



where 



F k - F k 



PL 
4 



E 



o 2 



N 

y — 

l= i^ k J^i dz ^ dz ™ 1 



(ek - ei) 2 



Nq- 



(25) 



(26) 



Note although F k , Li k , L 2k are still defined as in eqs.(18, 19) however the definition of Q is now slightly altered 
with fjv(g+i) replacing P^2 in cq.(20). Here fjv( 9 +i) is the joint probability density of q + 1 eigenfunctions, namely, 
Zi, Z 2 , .., Z q alongwith Zj (with j > q). Similar to previous case, the integral Q can again be approximated so as to 



' mn;k 



D- S (N - lY^x^mn - ELi z *mi z ni} r PNq- Here again x = 1 for n < (£ and 



express F k in terms of P/v g : Q r ri 
X ^ WC d ) for n > C d k . 

The above approximation for Q leaves the expression for F k in the same form as in eq.(23) however now hi = 
x(N — l)P Nq , and, h 2 = x(<5 mn — Ym=i z * ml z n i)PNq- The substitution of F k in F k gives the latter as a function of P^ q 
which in turn reduces eq.(26) in a closed form for Pjq q . The equation can then be used, by integrating over undesired 
components, to obtain the distributions of various combinations of eigenfunction components. 
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IV. DIFFUSION EQUATION FOR FLUCTUATION MEASURES OF EIGENFUNCTIONS 

The ensemble average of any measure of the eigenfunction correlations can be expressed in terms of P ( P oc P rq 
for a correlation function of r components of q eigenstates) . For example, the average of a measure, say C, describing 
the correlation among a set X of eigenfunction components can be written as 

/•OO 

(C(X;Y))= C(X-Y)P(X-Y)AX. (27) 
Jo 

where (.) denotes an averaging over various realizations of the sample. However the strength of the reproducible 
fluctuations of the correlations in different realizations of same complex system is of the order of the averages. As a 
consequence, a knowledge of just the averages is not enough and it is necessary to know the distributions of correlations. 

The ^-governed evolution of the distribution Pq of a measure C can be obtained by an integration of the undesired 
variables in eq.(14) (or eq.(18), eq.(25) as per requirement). As examples, we derive the evolution equations for few 
important measures in this section. The involved integrals are, however, quite tedious and analytical approximations 
seem necessary to reduce the equation in a closed form. As a check on our results, we study the Y — > oo limit of 
each measure. This limit corresponds to the flow of ensemble (1) (and ensemble (2)) to its steady state, that is, a 
Wigner-Dyson Ensemble. As a consequence, each measure is expected to evolve to its Wigner-Dyson limit as Y — > oo. 
We verify our results numerically too; the details are given in section VI. 



A. Distribution of Local Eigenfunction Intensity 

The distribution function of the local eigenfunction intensity i.e. the eigenfunction intensity u at a given basis state, 
say n can be defined as P u {u,e) — (X^fcLi — N\z n k\ 2 )5(e — e^)). The diffusion of P u as a function of A can be 
obtained from either eq.(14) with P oc Pin or eq.(18). For technical simplification, however, we choose the former 
and first study the evolution of the distribution Pn(x, e) of an eigenfunction component x = N 1 / 2 z n k = (u 1 / 2 ) at an 
energy e, defined as 

Pn(x,x*,e) = (6 f X) = J 5%P 1N {Z, E, Y) dE d?Z (28) 

whereof = 8(x— y/Nz n k)S^~ 1 (x* — y/Nz^ k ) and 5 e = S(e— e^) and dr = dEd^Z. The diffusion equation for Pi i(x, e) 
can be obtained by integrating eq.(14), with. P oc Piat; over the variables Gj and z n j^ j = 1 -> A, 



dPn = /3 2 
dY 4 

where 



d 2 d d(xG ) d(x*G ) 



2 + — — + 

dxdx* dx dx* 



+ L e Pn (29) 



G r (x,e)^ VJ J 8% J Zn f e [ P 1N dr, (30) 

with r = 0, 1 and / 5% S e [L E P 1N ] dE d^Z = L e Pn. 

Eq.(29) describes the sensitivity of the local intensity distribution to the energy scale e as well as various system 
parameters. As discussed in appendix A (see eq.(A7)), G r can be approximated as 



G r w »xo(N -\) l ~ T {N -\x\ 2 YPn{x)/til (31) 

with xo = M _1 f° r Ml x | 2 < 1 an( i Xo ~ \ x \ 2 for n\ x \ 2 > 1 where fj, — [e 2 f( Y ~ Y °'> — and Dk as the local mean level 
spacing at energy e^. A substitution of approximated G r in eq.(29) and an integration over e gives the energy-averaged 
local intensity distribution P x (x) = J Pu(x, e)de: 
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dP x 



4 



) d 2 [h 2 (x)P x ] , 5[/ii(x)P x ] , 0[/n(a:*)P x 



dxdx* 



O.r 



Ox 



with ft, 2 (^) = Xo( N ~ x 2 ), hi(x) = xo(N - l)x. Here A u = [iA with A = (Y - Y )/D 2 k . 
evolution of P x is governed by a rescaled parameter A u instead of Y. 



(32) 

Eq.(32) suggests that the 



<< N (thus A u = A), the above equation can easily be solved: P x (x, A\xq) oc e ^ x 1X °\ l 2 ^ x 7 ) 



For cases \x\ 

with 7 = q-< 3nk I 2 and Pr (xo) as the initial distribution. The steady state limit ^f- 
A — * oo. The solution in this limit corresponds to Wigner-Dyson case i.e. P x (x, A — * oo) oc ( 



of eq.(32) occurs at 
-0\x\ 2 /2 Qr ^ equivalently, 



Porter-Thomas distribution P u (u,A — > oo) oc u (' 3 - 2 )/ 2 e~ /3 "/ 2 p], @] (using u = |a;| 2 , which gives P„ = P.,.(2|x|) 

It is desirable to know the solution P x of eq.(32), or alternatively, P u for finite, non-zero A u and all ranges of u. In the 
diagonal representation of Ho, which corresponds to an initial distribution P Uo (uo, A = 0) = N~ 1 [S(u—l)+(N—l)S(u)], 
eq.(32) gives following short range behavior of P u : 



P u = [f5u/2f 2 - 1 



e -0u/2 _ 

! + -((/? + 2)//3 -{/3+ 2)V« + 0u/2) 



T((3/2) 
1 



exp 



(P/2)[-u+-u 2 + 



u<^ 2 



K- 1/2 < U < K- 1 



-2p na ( note re ~ ^ in i ar g e y_li m it and for D? - (/5iV) _1 ) 



(33) 
(34) 



where n = e 

The tail behavior of a distribution has a significant influence on its moments and the related physical properties. 
The asymptotic analysis of eq.(32) shows P u (u) to be a broad distribution: 



P u (u) 



exp 



-a u0 u 



1/2 



a un \vL n (nu) 



> 



(35) 



Here the coefficients are sensitive to system-specifics: a u o ~ 4qi/3 1 ( e ' 3JVA " — 1) ; a u \ ~ — AT/4, a„2 — (NP/16)e^ NAu , 
~ (-l) n (i/„/3 2 A^/4)e 2/3ArA " with z/„ decreasing as n increases. The decreasing coefficients alongwith alternate 
± signs lead to near-cancellation of higher order terms (with n > 2) in the exponent. Consequently, the tail is 
dominated by a log-normal behavior for the systems with large, finite A-strengths and a weaker than exponential 
decay in A — > limit. 

Eq.(35) indicates the existence of a log-normal asymptotic tail for the local eigcnfunction intensity of any complex 
system with finite, non-zero A u . A log-normal behavior of P(u) suggests a power-law behavior of its moments: 
(««) oc N~ dq [13 • Here d q is an effective dimension which can be different from a spatial dimension d. The form of 
P u (u) at finite A is therefore fixed by a spectrum of scaling exponents (as the moments can be used to recreate the 
distribution); the situation is termed as multifractal scaling. Further, as shown later, a log-normal tail of P(u) results 
in the similar behavior of the distributions of other related correlations and physical properties. Such a behavior 
has already been indicated for the physical properties e.g. conductance, density of states, local density of states and 
relaxation time etc. of disordered systems [3j. 

The significance of above P(u)-formulation is that here system dependence (other than size) enters only through 
one parameter, namely, A. This being valid for any complex system, modeled by eq.(l) (and eq.(2)), is thus applicable 
for disordered systems too. It is therefore relevant to compare our result with those obtained for disordered systems 
using other techniques (using renormalization group theory approach for dimension d — 2 + e, e < 1 [431 ]. and, by using 

Berezinski and Abrikosov-Ryzkhin techniques for strictly d = 1 cases [TJ IH3; the techniques predict a e~ aiU 
tail for d = 1, a log-normal tail for d = 2 and a log-cube tail for d — 3 case. However our technique predicts a 
log-polynomial behavior however dominated by log-normal term for all dimensions. 



B. Inverse Participation Ratio (IPR) 

The q th order inverse participation ratio I q of an eigenvector, say Uk, is defined as I q (k) — Ylj=i Wjk\ 2q - The 
physical meaning of I q can be illustrated by two limiting cases: (i) an eigenfunction with identical components 
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Ujk — N~ x / 2 corresponds to I q (k) — N x ~ q , and, (ii) an eigenfunction with only one non-zero component (say 
n th ) which gives Ujk — S n f. and I q (k) = 1. The case (i) corresponds to completely ergodic eigenfunctions covering 
randomly but uniformly the whole sample of volume V. The case (ii) corresponds to a wavefunction localized in 
the neighborhood of a single basis state. Thus I q , in general, is related to reciprocal of the number of components 
significantly different from zero and contains information about spread of a wavefunction in the basis space. For 
example, for a d-dimensional exponentially localized state, I2 ~ (a/(^) d , where a and £ are the lattice constant 
and localization length, respectively. Consequently, the typical value of I2 is a frequently used characteristic of the 
eigenfunction localization Q: = exp(ln/2) ~ N~ D2 with D2 a system dependent scaling exponent (also known 
as correlation dimension). 

The ensemble average of I q is related to g th moment of the distribution P x (x): (I q ) = N x ~ q J °° \x\ 2q P x (x)d^x = 
J IqPi dl q , The average inverse participation ratios can therefore provide information about the scaling exponents. 
As a consequence, it is useful to know the effect of changing system parameters on (I q ). Due to P (u) decay for the 
ranges [iu > 1, major contribution to (I q ) comes from the region [iu < 1. From eq.(32), it can be shown that 



^«ga<J,_i) -«*(/,). (36) 

where a — 2q + (3 — 2, t = N (3 + 2q — 2. Eq.(36) depends on two parameters, namely A and t which results in a 
different power law behavior for each (I q ), 



&(A)> 



a I (I q ^(r))e qtr dr 



(37) 



(I, 



9-1/ 



;Vi(A)) as 



For A — > 00, eq.(37) gives a correct steady state limit, namely, Wigner-Dyson behavior: (I q ) - 
^ 2 T q jN 1 ~ q for (3=1 and (I q ) = qlN 1 ^ 9 for (3 = 2. For finite nonzero A, (I q ) can be determined if 
well as some past information about the system (to choose it as an initial state which will give (I q (0))) is known. 
For example, for the systems where completely localized wavefunction dynamics is a valid physical possibility (e.g. 
disordered systems, mixed systems etc.), it can be chosen as the initial state which corresponds to (I q (0)) = 1; this 
gives (ii(A)) = 1, (I g (A)) w e ~^ NA for q < N. 

In general, the IPR fluctuations reflect the level to level variations of the spatial structure of eigenfunctions. In 
a complex system e.g. nano-system, however, the sample to sample fluctuations of the eigenfunctions also manifest 
themselves through IPR fluctuations which makes a knowledge of the /^-distribution over whole ensemble of samples 
relevant. The distribution Pj of I q of an eigenfunction, say Z^, with the components {^ n fe}fe=i,..iV is related to 
P oc P m : Pi q (I q ) = J <5i q PNi(Z k ,e k , Y) de k dr k with 5 Iq = 6(I q - £) n \z nk \ 2q ) and the volume element dr fc same 
as in eq.(20). The Y governed evolution of Pj can therefore be obtained from Eq.(18) for P oc Pni' 



DPj B 2 

-g^='-{X 1+ X 2 )+X, 



where X 3 = j 5i q [L E P N1 }dT k = and 

X x = J Si, [L lk + L* lk ]de k d' 3 Z k 
4 d 



1W hl5i aFl dr k , 



X 9 



2 / 5 Iq L 2k dT k 
8q 2 8 2 



(3 2 an 



SiF 2 dr k - 



4q(2q + f3-2) d 



f3 2 



ai n 



SiF 3 dr k 



(38) 

(39) 
(40) 
(41) 
(42) 
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with Fi 



Q 02 



nn;k 



)02* 
■* nn;k 



F 9 



Er, 



I ^mfc 



2(9-1) 



-^■n A; 



2(9-1)-* 



*n k ZmkQ% n ., k and F 3 = EnKkl^QiU where 



iife, ^2fe and Qmn-k are given by eqs.(19, 20). Using the approximate form (22) for Q™ n . k , F's can further be reduced: 

F 1 « q X (iV - 1)/,D 2 



X 

D 2 



£>2 



\ 2q - 2 -\z nk \ 2q ]P N1 . 



Nl 



(43) 



where x = 1 f° r A* < Cfe and x ~ A*/Cfe f° r M > C k - 

In general, the fluctuations of different moments (or measures) of the eigenfunction intensity need not be mutually 
independent. We can therefore define the joint distribution of two measures, say, h\(z), h 2 (z): 



P hl ,h 2 (hi,h 2 ) = J 5[h 1 -h 1 (z)]5[h 2 -h 2 {z)]P N1 de k d f3 Z k 



(44) 



The above definition along with the equality £„ \z nk \ 2 ^~^ = £ m n \z m k\ 2q \z nk \ 2 ^ - E m ,n;m?n \z m k\ 2q \z n k\ 2( " ! 
gives 



2(2g-l) 



P N1 de k d f3 Z k 



I 



/gVlP/^^VOdVi- / WP /9 , W (/„W)dW 



^ 9 - K vpp h 



(45) 



where is a measure of the correlation between the intensities localized at two different basis sates: W q — 
Em n-m^n I z mk \ 2q \ z nk \ 2 ^ 9 ~ ^ ■ The 2 nd equality in eq.(45) is obtained from first by replacing l q -\ and W by their 
typical values; the superscript typ over a variable R indicates its typical value: R typ = exp(lni?).) Using eqs.(44,45), 
the terms X\ and X 2 can be rewritten as the functions of I q and Pj q which in turn leads to 



dPi a 



dA ip ' ~* dig 2 



with a, t same as in eq.(36) and 



K lp = q X A. 



(46) 



(47) 



Note, the above equation alongwith the definition (I q ) = J I q Pi dl q again leads to eq.(36). 

The behavior of Pj q in different I q regimes can now be probed by analyzing eq.(46), using completely localized 
eigenstates as the initial state. The behavior varies from an exponential decay for small- I q regime to log-power law 
decay for asymptotic tail regime of I q : 



P Iq ~ exp[-a i0 /g] I q <e- X l^ 



~ I. x exp 



M 



-E^ln"(/ g //* w ) 



n=l 



(48) 
(49) 



with e w 2.72, a i0 w t(l - e"*^)" 1 /2q and a in - (-l)"2(/* w 1 ) n /(2 g /*2'P)™ (valid for q < N) and M as a large 
integer. Note the alternate ± signs of terms with increasing powers lead to convergence of the series in the exponent. 
However the tail is dominated by increasingly higher powers of the logarithmic term as I q increases above its typical 
value. For example, for e~ 1 I t q yp < I q < I q yp , n = 1 dominates the exponent and Pi q behaves as a power-law. Similarly, 
the tail shows a log- normal decay for regime I* vp < I q < eIq VP . 
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Eq.(46) depends on more than one parameter, namely, A ip as well as size-dependent parameters (appearing through 
t). This suggests an absence of single parameter scaling in IPR distributions. However, as suggested by eqs.(48, 49), 
it seems possible to define a single parameter locally (that is, different single parameters governing different IPR 
regimes). Further note that the asymptotic behavior of Pi q is sensitive to A-strength and is therefore system-specific. 
This result also agrees with the NLSM-result obtained for disordered systems [13]. 



C. Pair-Function w(r,r') 



The measure contains important information about the spatial correlations between components of an eigenfunction 
Zj at two different basis points of the sample and at an energy e: w(n,m) — \z n jZ m j\ 2 ; (equivalently w(r,r') = 
\zj(r)zj(r')\ 2 in a continuous basis e.g. coordinate space r). In the localized phase, the asymptotic behavior of 
lnui(r, r') at |r — r')| — > oo determines the rate of exponential decay of the eigenfunction amplitude. It is also useful 
for many physical applications e.g. in determination of the form factor of resonance scattering in the c omp lex nuclei 
or the resonance conductance of the quantum dot with point contacts in the coulomb blockade regime [15| . 

The distribution P w . e = (X^ 6(w — \z n jZ m j\ 2 )5(e — ej)) of the correlation between n th and m th components of an 
eigenfunction, at an energy e, is related to Pnx- P w (w, e) = ^ J 5 w 5 e jPNi(Zj, ej, Y) drj with S w = 6(w— \ z n jZ m j\ 2 ) 
and drj same as in eq.(20). Consequently, its rate of change with respect to Y can be determined by eq.(18), 



BP B 2 

^=NL e P w + ^(A 1+ A 2 ) (50) 

where 

^ = | E ^ w J S ™ 6 ^ [Qnn-J + dr, (51) 

A2 = 4 E ^2 w [ Med [Fx + F 2 ] de j d^Z j 5 ™ 6 ^ ^ + 2 ^ _1 ^] dr * ^ 

3 3 

with Fx = \z mj \ 2 Q 12 n , 3 + knj| 2 Q^ m;J and F 2 = z* j z m jQ™ n . j + ZnjZ^Q^.j. Eq.(50) is derived by first using eq.(18), 
followed by repeated partial integration. Note J S w 6 e j[LEPNx]dejdZj = NL e P w . Within approximation (20) for Qs, 
Ax, A 2 can further be simplified which on substitution in eq.(50) give the diffusion of P w in a closed form: 

BP 8 2 B 

7^ = «-2 N n i - 4u; )l P ™-7T [ fi 2 - H + NL e P w (53) 
oA w aw A aw 

where fix = \z m] \ 2 yp + \z n A 2 yp + 2\z nj \ 2 yp S nm = 2(1 + S nm )uf p , with u t J vp = exp[(\aUj)] as the typical local intensity of 
the j th eigenfunction, ft 2 = {P/2)[\z mj \ 2 typ + \z nj \ 2 yp + {4/P)\z nj \ 2 yp 5 nm \ = {(3 + 28 nm )uf p , A w = 2 X A and 6 = (3N + 2. 
The last term on the right of eq.(53) can be removed by an integration over energy e, leaving us with an evolution 
equation for energy-averaged P w . Note the energy averaging of eq.(53) for case n — m corresponds to eq.(46) for 
P(J 2 ) (as = -&). 

Exploiting the similarity of the form of energy-averaged eq. (53) to eq. (46) , the behavior of P w (w) in different regimes 
can again be given by eqs. (48,49) after following replacements (everywhere in the equations): I q — > w, (Xi n — > a wn 
where a wn ~ (-l) n 2 1 ~ 2n fl%(w typ )~ n for n> 1 and a w0 « 6(1 - e - ^) -1 / 4 - Thus P w (w) decays exponentially for 
small w ranges (w < w typ /e): P w ~ exp [— a w ow\. It shows a power-law behavior for regimes e~ 1 w typ ^ w > w typ 
P w ~ w - 1 e ~ a ™i in ( KW ) 5 a log-normal decay for regimes w typ > w > ew typ . Such a behavior was predicted by non-linear 
sigma model studies of quasi 1-d disordered wires too [H . 

Eq.(53) can be used to study the behavior of various moments of the distribution of pair-correlation. For example, 
for average behavior of w, that is, (w) — J wP w (w; A) duo, eq.(53) gives its A evolution. The evolution equation turns 
out to be of the same form as eq.(36) (with q = 2) with following replacements: (I q ) — > (w), a — > fl 2 and 2A — > A w 
(note (I q -x) = 1, t-> b for q = 2). It can be solved to show that (uu(A w )) = e~ bA ™ [(w(0)) + (n 2 /b)(e bA ™ - 1)]. A 
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choice of (w(0)) = S m i (corresponding to localized regime) gives (w(A w )} = <5 TO ;e ™ + (fi 2 /i>)(l 
which is analogous to the result obtained for disordered systems (by NLSM techniques); see 



-6A„ 



{\-k)/N 



D. Correlation Between Eigenfunctions At Two Different Energies 



Critical point studies of many systems indicate the presence of multifractal structures among eigenfunctions. The 
multifractality suggests that the wavefunction is effectively located in a vanishingly small fraction of the system 
volume. However such extremely sparse wavefunctions can exhibit strong correlations if they belong to neighboring 
energy levels; the correlations therefore preserve the level- repulsion despite the sparsity of the wavefunction. Thus, for a 
complete analysis of level-statistics and associated physical properties, a knowledge of correlations among eigenfunction 
is very important. The correlations are also used in the analysis of many other physical properties e.g. for the 
measurement of the linear response of the system, or, to determine the fluctuations of matrix elements of some 
operator in a given basis. This information is useful in studies of the effect of a particular interaction on the statistical 
properties of the system e.g. effect of electron-electron interaction on a single particle disordered system. 

The correlations between components of two eigenfunctions at different energies can be described as cr(n, m, e k , e;) = 
\ z nkZ m i\ 2 ) (equivalently, in a continuous basis: a(r,r' ,e,e') = \ip e (r)ip e > (r')\ 2 ). The distribution P a of the correla- 



2 between n th and m th components of the eigenfunctions Zf. and Z[, respectively, is related to 
= Lk,i/ ^<T <5 e,k^e+ li ;,i P N2(Zk, Zi,e k ,e 1; Y) UIj=k,l dr j) with S <j = S ( a ~ \ z nk z m i\ 2 ), d.Tj same as in 
\e k — ei\ as the energy difference between two states. Using eq.(25) and proceeding as in the case 
of P w , the diffusion of P a (cr) can be shown to be described by the equation 



tion a = \z nk z nl 
cq.(20), and, u> 



dPa_ 



}2 r 



8a 2 



a(Qi — a) 



Oct 



ba) 



L e + L 



Prr 



(54) 



typ 



\Z r k\ 



typ 



2\ z rk\ 2 t y p 8kl8rs), ^2 



Z rk\ 2 tvp +\z 2 Hyp ) +Mz rk \ 



typ 



SkiS rs ]/2 and 



where A a = 2^A, tii = (\z s 

b = f3N + 2). Note eq.(53) is a special case of the above equation (as P(w) = P{cr) for k = I). 

The energy averaging of eq.(54) once again leads to an equation similar in form to eq.(46). Exploiting the analogy, 
we again get three different regimes for Pcr(a): (i) ~ exp [— a a Q(j] (for a < e~ 1 (T typ ), (ii) ~ cr _1 exp[— a<ji\n(a / a typ )] 
(for e- x a tw >a> a typ and (hi) ~ o-^xp [-a^ 2 {a / a^)} for ea^a > a typ where a<? m 6(1 - c- bA -)-V4 and 
a an - (-l) n 2 1 - 2n f2J(cr t w)- n for n > 1. 

The A dependence of average behavior of (a) can now be derived by multiplying eq.(54) by a and then integrating 

a(0)) + ({l 2 /b)(c bA ° - 1) . 
= which gives (cr(A)) w 



over a; the equation again turn out to be same as the q = 2 case of eq.(36) after following replacements 
a — > ^2 and 2A — > A^. Solving the so-obtained evolution equation gives (a(A,j)) = e~ bA " i 

The choice of a localized initial state (e.g. an insulator at A a = 0) corresponds to (<x(0)) 





N 



In this paper, we have considered only two point correlations. The other correlations e.g. (z* k z r iz s kz* sl ) related 
to linear response of the system or, higher order ones e.g. (|2 r fc| 4 |z s fe| 4 ) related to IPR fluctuations, can similarly be 
determined using eq.(25). 



E. Local Density of States p(e,j) 



The local density of states or the spectral function, defined as p(e,j) — Y] \Uj n )\ 2 8{e — e n ), is an important measure 
of localization. This is because it counts the eigenstates U n having appreciable overlap with (or equivalently, located 
close to) the site j. Note this is distinct from the global density of states p(e) which counts all the eigenstates at the 
energy e irrespective of their location in space. The measure p(e,j) is of special interest as it gives information about 
the decay of a specific unperturbed state into other states due to interaction. The width of the LDOS defines the 
effective "life-time" of the unperturbed basis state. Its distribution is an experimentally accessible quantity related 
to the position and form of NMR line [|[ . 
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The probability density P p {p) of p(e,j) is related to Pin- P p (p) = J <5pPiN(E, Z, Y) dE d^Z where S p — S(p ~ 
|-Zjn| 2 <5(e — e„)) (note here Z = {zj n } n —i ..jv). The diffusion of P p due to changing system parameters can again 
be studied with the help of eq.(14) for Puv- 

OP 

= L E P P + B + B* (55) 

with B = B* = j dp L Z P 1N d^ZdE = ^ + £§*z. Here the 2 nd form of B is obtained from the I s * by a substitution 
of three terms of L z (eq.(15)) in the integral, and, a subsequent partial integration which gives 

Fl . Js f 

F, . / { „ 

As in the case of the integrals Q and G (see Appendix A), the dominant contribution to the integrals F\ and F2 
comes from the regions where the exponent term in Pxn, that is / = (p/2) J2 m<n |e„ — e m \ 2 \Uj n \ 2 \Uj m \ 2 < 1. This 
can occur under two conditions: 

(1) p(a) < 1: Here er = |L^- n | 2 |£7j m | 2 describes the correlation between two different eigenfunction components in 
the same basis state. Under the condition, a neighborhood of the order of mean level spacing can contribute to the 
integral over e-variables i.e \e n — e m | 2 ~ D 2 , 

(2) p(a) > 1: In this case, / 3> 1 only for those regions where |e„ — e m | 2 ~ P/ 2 / (p(a)). Note however in both cases, 
almost entire eigenfunction space can contribute to the integral. 

Thus Fi and F2 can be approximated as F\ « xF>~ 2 [p(l — p)-Pjo] an d F2 w NxD~ 2 [p — (p)] P p where \ = 1 if 
/i(cr) < 1 and x ~ A t ( <7 ) f° r / i ( (T ) ~ 1- The approximate forms of Pi and F2 can now be used to rewrite B as a function 
of P p which on substitution in eq.(55) gives the diffusion equation for P p : 



\ Z jn\ 2 \Zj m \' 



J 



<J(e - e n ) (1 - (5(e - e m )) 



Piat d^ZdF 



(56) 



E 



(e r , 



■5(e - e n ) 



(57) 



m; ^ ieP " + ^ w ^ )] p ' + [ " - ( " )] p ' (58) 

where A p = %A. Note that the above equation is analogous in form as the eqs.(46, 53, 54) of Pj , P w and P CT , 
respectively. This similarity is reflected in both short as well as long-range behavior of P p : (i) ~ exp [— a p op] (for 
P ^ P typ ), (h) — p _1 exp— a p i\n(Kp) (for n^ 1 > p > /9 tl,p ), and, (iii) ~ p _1 exp [~a 3 ln 2 (/«/?)] for (<r > ft -1 ) where 
a P o « 6(1 - e- 6A ")-V4 and a p „ - (_l)n 2 i-2n( (0 t OT )-n_ 



V. THE PARAMETER A 



The set of Eqs. (14, 18,25) provides a common mathematical formulation for the eigenfunction-statistics of various 
complex systems modeled by eqs. (1,2); here the information about the system enters only through Y. As shown 
explicitly in [4!|, the same Y also enters in the common mathematical formulation of the eigenvalue-statistics of 
ensemble (1) (and (2), see (EH); this is implied by eqs. (14, 18, 25) too. However, as discussed in 0, [5l|, the evolution 
of the n th order eigenvalue correlations (n > 1) as a function of Y, is abrupt in large N- limit; a smooth crossover can 
only be seen in terms of a rescaled parameter A e where 



A e (e,y)=A=I^° (59) 

with D^(e, Y) = D((/L) d as the local mean level spacing, D(e, Y) as the mean level spacing of the full spectrum and C, 
as the correlation/localization length for a d-dimensional system of length L (N — L d ), at an energy e and parameter 
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Y (with Y as its initial value). Thus A e for various systems e.g disordered systems, mixed systems, systems with chiral 
or particle-hole symmetry etc. can be calculated by a prior knowledge of system parameters; (e.g. see for the 
calculation for Anderson and Brownian ensembles). As A e increases from zero to infinity, the level-statistics changes 
from its initial state (with Y = Yo) to that of Wigner-Dyson ensemble. For example, let the initial state corresponds 
to insulator limit of disordered systems or integrable limit of mixed systems; both limits show Poisson level-statistics 
0, S S3, [EH, [53)]. A variation of system parameters changes A from zero, causing diffusion of levels towards Wigner- 
Dyson steady state. According to A formulation, the level-statistics, for the system parameters resulting in finite A, 
is then an intermediate point of Poisson — > Wigner-Dyson transition. The prediction is in agreement with previous 
works on the two systems 0, HH [EH [El, [E3] ; (note Wigner-Dyson statistics corresponds to metal and chaotic limits 
of disordered and mixed systems, respectively [13, IHj]). 

For later reference, it is worth reviewing the role of A e , that is, A in locating the critical point of level statistics. 
As both \Y — Yq\ as well as the local mean level density are functions of N, the latter can affect A significantly. As a 
consequence, the size N plays a crucial role in determining the level statistics in the critical regime. For finite systems, 
the eigenvalue statistics smoothly approaches one of the two end points, namely, A — > or A — > 00, with increasing 
system size. The variation of A in infinite systems, however, may lead to an abrupt transition of the statistics, with 
its critical point given by the condition A = size independent (see ref.[l7j for the definition of a critical distribution). 
The finite, non-zero A strength, say Acriticai, at the critical point, results in an eigenvalue-statistics different from 
the two end points. Note, however, that the existence of a critical point or its absence depends on the relative size- 
dependencies of \Y — Yq\ and the local mean level spacing. If the size-dependence of remains different from that 
of \Y — Yq\ under all complexity conditions, A will never achieve a finite non-zero value in infinite size limit. As a 
consequence, such a system will not show a critical behavior of eigenvalue-statistics. For example, as discussed in (Hoj 
for a d-dimensional Anderson Hamiltonian of linear size L, A turns out to be size-independent only for d > 2. The 
A-formulation, therefore, indicates the lack of metal-insulator transition for dimensions d < 2 which is in agreement 
with several studies of previous years. 

The connection of the eigenvalue fluctuations to those of eigenfunctions suggested A as the evolution parameter 
for the eigenfunctions correlations (of order n > 1) too. As shown in section III, the evolution parameters A measure 
of various eigenfunction fluctuation measures are indeed functions of A: A measure = /(A). Here /(A) oc A on short 
length scales and /(A) ~ Ae~ aA in the tail regime. 

The parameter A, being a function of the distribution parameters of the matrix elements, is sensitive to changes in 
the system parameters; this is due to latter's influence on the uncertainties associated with system-interactions. Some 
examples of such system parameters are disorder, dimensionality, boundary and topological conditions, system size etc. 
For example, the presence of disorder randomizes the interactions, with degree of disorder affecting the distribution 
parameters h, b and consequently A. The dependence of A on the dimensionality and boundary conditions originates 
from their influence on the basis connectivity i.e degree of sparsity of the matrix which is reflected in the distribution 
parameters h, b. For example, for nearest neighbor hopping and hard wall boundary conditions in (i-dimensions, the 
matrix element Hjk ^ only if j = \k ± L d_1 | (with L as linear size). The variance hjk of the distribution p(Hjk) is 
therefore finite only for j = k or \ j = k ± L d ~ 1 \ and is zero for all other j, k. The information about dimensionality in 
A also enters through the local mean level spacing which depends on the correlation volume £ d . (See also where 
the dependence of A on system parameters is explained by considering an example of Anderson Hamiltonian) . 

The system size N is another important parameter which affects the evolution of the measures significantly. As 
shown in section IV, it appears independently as well as through A measure in the evolution equations which suggests 
a two parametric dependence, namely, A measure and N (separately) of these measures. As a consequence, even at the 
critical point of level statistics, the eigenfunction statistics remains sensitive to size N. This in turn results in a mul- 
tifractal behavior of the eigenfunctions at the critical point of any complex system, modeled by eqs.(l,2). The scaling 
exponents at the critical point, referred as critical exponents or multifractal dimensions, depend on system parameters. 
In finite size systems, changing system parameters can change A (and therefore A measure ) continuously between 
and 00 which may lead to intermediate stages of varying degree of multifractality. However, the physically interesting 
cases usually correspond to infinite sizes where A takes only three possible values, namely, A = 0, 00, A cr iticai] for 
these cases therefore only one multifractal stage, that is at the critical point corresponding to A cr iucai, can exist. As 
Acriticai is sensitive to system-specifics, the critical (multifractal) exponents can vary from system to system. Note, 
as already mentioned above, the occurrence of critical point and, therefore, a multifractal behavior of eigenstates is 
not a necessary feature of all infinite size complex systems. 

The A-governed diffusion equations, derived in section IV, are valid for arbitrary initial conditions at A = (which 
implies A x = 0) and their solutions P x (X,A x \X o ,0) describe the probability of the measure, say V, at A x for a given 
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initial state of X = Xq. Thus P is subjected to an initial constraint limA — > QP(X, A X \X , 0) = 5(X — X ). By 
integration of the solution over the distribution of initial values Pq(X o ,0), one can recover P(X,A X ), that is, the 
distribution of measure X for a system with complexity parameter strength A: 



The eq.(60) implies that the statistics evolved in "time" A x is sensitive to the collective behavior of system pa- 
rameters contributing to A x and the initial distribution only. The latter can always be chosen same for the systems 
operating in the matrix spaces of similar type e.g. same symmetry conditions; (the initial values of their Y parame- 
ters need not be equal). Thus if both A,B operate in the matrix spaces of same type, their behavior at the system 
parameter strengths which lead to A x ,a = A x ,b = t would also be same (although they may show different behavior 
between < A x ^a, A^.b < t). This implies a great deal of universality among systems of widely different origins of 
complexity. For example, consider the cases of a three dimensional disordered system, say A, and a clean, closed 
quantum dot, say B. In the first case, A x _a = A x _disorder is a function of disorder, hopping strength, dimensionality, 
boundary condition etc. In the case of a dot, A Xj b — A x dot is a function of shape as well as size. It is well-known that, 
in strong disorder limit and for circular shape, respectively, both systems show localized wavefunction dynamics and 
same statistical behavior of the eigenfunctions and eigenvalues of the Hamiltonians. Reducing the degree of disorder 
or change of shape of the dot from circle to stadium type results in a transition from localized to delocalized dynamics 
of the wavefunctions. The statistics in the intermediate stages during the transition for each case is governed by 
the respective A x strengths. If, however, the A Xt dot = A x> disorder at some shape parameter and disorder strength, 
respectively, our analysis predicts a same statistical behavior for both systems. The implication can also be extended 
to classical systems e.g. stock market fluctuations which are analyzed by statistical studies of the correlation matrix 
of stocks (Note, the correlation matrices of classical systems are, in general, non-Hermitian; however, as shown 
in [4fj| . the A formulation remains valid for non-Hermitian version of eq.(l)). Here a very weak interaction among 
certain stocks due to various socio-economic conditions results in a localized dynamics of the eigenfunctions. The 
changing conditions may lead to a more homogenized interaction of some of the stocks, thus introducing a transition 
from localized to delocalized wave dynamics. In this case, A is a function of the socio-economic parameters (SEP). 
However if A XtStoc k = A x ^dot for some combinations of SEP and dot parameters, respectively, the spectral and strength 
fluctuations in correlation matrix of the stock market and the Hamiltonian of quantum dot will show same behavior. 
Note the analogy of statistical behavior of the eigenvalues and eigenfunctions among the three systems, mentioned 
above, has already been numerically verified in delocalized waves limit A oo dUGl. . 

The above universality makes A formulation useful as it can be exploited to obtain the statistics of a complex system 
if the same information is available about another system (under same symmetry conditions) by another method. 
For example, for Anderson type disordered Hamiltonian, the distributions of many measures are known by non-linear 
sigma model techniques. The formulations can then be used for the complex systems e.g stock markets undergoing a 
transition from localized — * delocalized wave dynamics; one just needs to replace A(Anderson) by that of the system. 

The formulation can also be used to search for the system conditions leading to a critical state or multifractal 
wavefunctions of various complex systems. For example, the suggested modelling of mixed systems by eq.(l) would 
imply the possible existence of a critical point of level-statistics in the systems and multifractal eigenstates. The 
intuition suggests that the occurrence of such a point may correspond to breaking of the last KAM curve, thus allowing 
classical diffusion or derealization of the dynamics above the critical point and localization below it; however it needs 
to be further explored. The critical A can then given by the critical value of system parameter leading to last KAM 
curve breaking. 



For numerical analysis, we choose three different ensembles (for both cases j3 — 1,2); the choice is dictate d by 
the reason (i) the ensembles are prototype models of many physical systems related to different areas [3j, |g, |{§ |35j , 
(ii) a comparative study of the eigenvalue fluctuations of these systems has already been carried out, with their A 
parameters and other results given in [49| : 

(i) Critical Anderson ensemble (AE): 

(a) Time-Reversal case AE t : we analyze cubic (d — 3) Anderson lattice of linear size L (N = L d ) with a Gaussian 
site disorder (of variance W 2 /12, W = 4.05 and mean zero), same for each site, an isotropic random hopping between 




(60) 
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nearest neighbors with hard wall boundary conditions [3j, 1 501 ] . The ensemble density in this case can be described by 
eq.(l) with hkk = W 2 /12,h k i = f(kl)/12,bki = where f(kl) = 1 for {kj} pairs representing hopping, f(kl) — » 
for all {k,l} values corresponding to disconnected sites. A substitution of above values in eq.(8) gives Y which 
subsequently gives A by eq.(59): (see eq.(19) of [49|). 

K(E,Y) = la-aolF 2 ^"^" 1 (61) 

with a — ao = 1.36 and F(E) = 0.26e~ B / 5 (see section V of 03)- (Note, for later reference, F(E) is the mean level 
density: F(E) = (ND)' 1 ). 

(b) Broken Time-Reversal case AE rat : we analyze cubic (d = 3) Anderson lattice of linear size L (N = L d ) with a 
Gaussian site disorder (of variance W 2 /12, W = 21.3 and mean zero), same for each site, an isotropic non-random 
hopping t = 1 between nearest neighbors with periodic boundary conditions [3j, 1 501 ] . The time-reversal symmetry 
is broken by applying an Aharnov Bohm flux <fi which gives rise to a nearest neighbor hopping Hki — exp(icf>) for 
all k,l values related to the nearest- neighbor pairs [58| . The flux is chosen to be non-random in nature, that is, 
(cos 2 (0)) = W\ — 0, < sin 2 (0) >= Wi = and (cos(<^>)) = t\ = 1, < sin(0) >= i 2 = 1. The ensemble density in this 
case can again be described by eq.(l) with hkk = W 2 /12, bkk — 0, hki- s =W S = 0, bki- s — t s f{kl; s) where f(kl; s) = 1 
for {k, 1} pairs representing hopping, f(kl; s) — > for all {fc, 1} values corresponding to disconnected sites. The A for 
this case is still given by eq.(61) (except for a factor /3 _1 ); however now a — a = 5.43, F(E) = 0.016e~ £ Z 400 (see 
section V of (U). 

(ii) Critical Power Law Random Banded Ensemble (critical PRBM or PE): As mentioned in section II, PRBM 
ensemble was introduced as a possible model for the level statistics of Anderson Hamiltonian [34]. It is defined 
as the ensemble of random Hermitian matrices with matrix elements as independently distributed Gaussian 
variables with zero mean i.e < Hij >= and a power-law decay of the variances away from the diagonal (3. fiol. l34j : 
< \Hij- s \ 2 >= a(\i — j\) with function a(r) decaying as r increases. 

The PRBM ensemble with specific choice < \H lj;s \ 2 >= G^ 1 [l + (\i - j\/b) 2 ] ~\ Gy = (3(2 - S v ■.) and Gy = 1/2 
(referred as critical PRBM or PE in this paper) leads to a critical behavior of eigenfunction and eigenvalue statistics 
at arbitrary values of the parameter b and is believed to show all the key features of the Anderson critical point, 
including multifractality of eigenfunctions and the fractional spectral compressibility [H, HH . The ensemble density 

in this case corresponds to eq.(l) with bki = 0, and, hu-s = G^ 1 [l + (\k — l\/p) 2 ] . The corresponding A can be 
shown to be given by (see section VI of |49j]). 

A p (p, E) = a ] ; 1 f(p)F 2 (E)C 2 N- 1 . (62) 

where a p = 2N(N + 2-/3), f(p) = J2r=i( N - + (p/ r ) 2 \- 

(iii) Critical Brownian Ensemble (BE): A Brownian ensemble can be described as a non-stationary state of the matrix 
elements undergoing a cross-over due to a random perturbation of a stationary ensemble by another one [1, [II], [H, . 
For example, in the case of Hermitian operators, a Brownian ensemble H can be given as H = \J~] (Ho + ^V) (with 
/ = (1 — A 2 ) -1 ); here V is a random perturbation of strength A, taken from a stationary ensemble [33] e.g. Wigner- 
Dyson ensemble, and applied to an initial stationary state Hq (see also [II]]). Here we consider a specific class of 
BEs, namely, those appearing during a transition from Poisson — » Wigner-Dyson ensemble, caused by a perturbation 
of the former by the latter (that is, taking Hq, V as Poisson and Wigner-Dyson ensemble respectively). As, in 
above two cases, this transition also results in a change of localized eigenstates to delocalized ones. The BEs related 
to the Poisson — > Wigner-Dyson transition can be described by a N x N ensemble H represented by eq.(l) with 
mean (H M ) = b ki = 0, the variance (Hl k . s ) = h kk[S = (2-f)^ 1 and (ff 2 ; . s ) = hu-s = [47(1 + p)}^ 1 for k ^ I. with 
(1 + fi) = (A 2 /) -1 ; here H = Hq for A — > or /i — > oo. As mentioned in section II, the ensemble density in this 
case has the same form as for Rosenzweig-Porter (RP) ensemble [33|; it can also describe an ensemble of Anderson 
Hamiltonians with very long range, isotropic, random hopping. Further, as discussed in [Ii| , the special case /i = cN 2 
corresponds to the critical BEs; their mean level density is given as F(E) = (ir)~ 1 / 2 e~ E and 

A b (E) = (l/4c7r7)e~ £2 . (63) 

Our aim is to show that the behavior of an eigenfunction fluctuation measure of AE, BE and PE is analogous 
at system parameters which lead to a same A measure value for all the three cases. Using the latter as a condition, 
we can obtain the desired system parameters in each case (that is, p for PE and c for BE for a given AE). As A 
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for the three cases is energy dependent, the fluctuation measures should be compared at precisely a given value of 
energy. For numerical analysis, however, one needs to consider averages over an energy range AE which should be 
sufficiently large in order to improve the statistics. On the other hand, choice of a very large AE will lead to mixing 
of different statistics (in a range A A oc 5E). As a consequence, one needs to consider an optimized range of AE. In 
our simulations, we analyze large ensembles of about 1400 matrices of size N — 2197. We choose AE to be about 
10% of the bandwidth, at the band center E = which gives approximately 3 x 10 5 levels for each case. As the chosen 
AE corresponds to a 1 % variation of the density of states, it avoids mixing of different statistics. 

As discussed in previous section, the eigenfunction fluctuations are influenced by both A as well as system size N. 
To compare A measure dependence of a fluctuation measure (of eigenfunctions), therefore, same system size should be 
taken for all systems under consideration. As examples, here we consider distributions of three measures, namely, 
local eigenfunction intensity P u (u), inverse participation ratio Pjfa) and pair correlation function P w (w) for the three 
systems under time-reversal symmetry (j3 = 1) i.e . AE t , BE t , PE t . As, for P u , 

_ f ,(Y-Y X 2d / F V 

yV « i\T9. n9 ~ i-tvv ' V 04 J 



N 2 D 2 \I* VP . 

(with C, d ~ (l2 VP )~ 1 )i the BE and PE analogs for the intensity distribution of AE t can be obtained by the condition 
1*2 a I Fa — I^b/Fb = I^p / ' F p . This requires a prior information about ll vp . Our numerical study for various sizes 
of the three systems shows that, for each case, w IN~ D2 with / and system dependent. The numerical 
information about and F can now be used to obtain the parameters p and c for PE and BE analogs of AE t for 
P u case (i.e PE and BE with the ratio /F same as for AE); we find p — 0.4, c = 0.02. The figure 1 shows the 
distribution P u (u'), u' — [lmt — (lnu)]/(ln 2 w), for the BE t case (c = 0.02) and PE t case (p = 0.4) along with AE t 
case; the close agreement among the three cases confirms our theoretical prediction. This is also confirmed by the 
comparison of Piling) and P w (hnv) for the three systems, shown in figure 3 and figure 5, respectively. Here again 
the parameters p and c for PE and BE analogs for both measures are obtained by the relation A/ a = Ajj, = A/. p 
(similarly for w). 

The above numerical analysis is repeated also for the case of AE in a magnetic field and its BE and PE analog; the 
results for the three measures, shown in figures 2,4 and 6, further support our claim: the eigenfunction fluctuations 
of different complex systems show same behavior if their complexity parameters and sizes are equal. It is worth 
reminding that behavior of the eigenvalue fluctuations is governed only by the related complexity parameter ( that is, 
no independent influence of size). The details of analytical and numerical evidence about the eigenvalue statistics are 
already published in [49l. l50l. l5l| . However for the sake of completeness and to convince the reader, we include here 
the numerical analysis of an eigenvalue fluctuation measure, namely, nearest neighbor spacing distribution P{S) for 
the three system (for both f3 = 1,2 cases) at parametric values leading to A e a = A e ^ = A EyC (where A e = A); the 
plots shown in figures 7,8 reconfirm the claim about eigenvalue statistics. 



VII. CONCLUSION 



In the end, we summarize our main results. Our analysis of the eigenfunction correlations of complex systems 
indicates a two parameter dependence, namely, complexity parameter A and system size N, of the distributions of 
eigenfunction components. The independent appearance of size parameter (besides through A) seems to suggest a 
lack of finite size scaling in eigenfunction distributions and an absence of their critical limit. This is in contrast 
with the behavior of eigenvalue distribution which shows a single parametric scaling as well as a critical limit if the 
condition limA^ — > oo A = finite is satisfied by the system. Note the above implies the size-dependence of eigenfunction 
correlations at the critical point of level-statistics too. 

We have also studied the distribution of a few important measures of eigenfunction correlations e.g. local density 
of states, pair correlation function etc. We find that the form of complexity parameter governing an eigenfunction 
fluctuation measure is sensitive to its nature (e.g. A u for local intensity distribution, A/ for inverse participation 
ratio distribution etc). This is again different from the eigenvalue fluctuations (except for level density) which are 
all governed by the same complexity parameter, namely, A e = A. Our analysis indicates a log- normal behavior of 
the asymptotic tails of the distributions at finite A-strength. In context of disordered systems, a similar behavior 
was predicted by other studies using different techniques e.g. Berezinski and Abrikosov-Ryzkhin technique (for 
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one dimension) and by non-linear sigma model (for higher dimensions) Q. However the complexity parameter 
formulation suggests the existence of such a tail-behavior and multifractal eigenfunctions for almost any complex 
system, irrespective of the origin of complexity, if the parameter N A measure is finite. A recent numerical study of the 
eigenfunction of the correlation matrix of stock prices confirms the suggestion in case of stock market [5|, [28[ . As finite 
A corresponds to the critical point condition in infinite size systems, a log-normal tail-behavior seems to be associated 
with the existence of a critical point (and vice- verse). The above study can thus be used to search and predict the 
critical stages of other complex systems e.g. stock market, brain etc. 

In this paper, we have considered the cases modeled by generalized Gaussian ensembles with uncorrelated matrix 
elements as well as a wide range of non-Gaussian ensembles with correlated matrix elements. The latter are suitable 
models, for example, for disordered systems with varying degree of particle-particle interactions. In context of dis- 
ordered systems, therefore, we expect a same statistical behavior of a measure for both the cases, namely, with or 
without particle interactions, if the strengths of their parameters A measure are equal. This suggests the sensitivity of 
the statistical behavior of a disordered system to degree of its complexity only (measured by complexity parameter), 
irrespective of the origin. The statement is expected to be valid for correlated and uncorrelated cases of other complex 
systems too. This in turn would indicate the existence of an infinite family of universality classes, parametrized by 
A, of statistical behavior among complex systems. 



APPENDIX A: CALCULATION OF INTEGRALS Q™ n . k AND G r 



The integral Q^ n . k defined by eq.(20) can be rewritten in terms of p(H) as 



Q™* = £ / h P(H, Y) dH (Al) 

To express Q in terms of Pjvi, it is necessary to write p(H) in eigenvalue-eigenvector space i.e. {A, U} space. 
The steps can briefly be given as follows. The solution of eq.(6) for arbitrary initial condition, say Hq at Y = Yq 
can be given as p(H,Y\H ,Y ) cx exp[-(a/2)Tr(iJ - t]Hq) 2 } with a = 7(1 - rj 2 )^ 1 and 77 = e' lY . Without loss of 
generality, the basis space for H can be chosen as the eigenvector space of Hq; The initial ensemble Hq in this basis 

consists of diagonal matrices. For simplification, consider the initial distribution given by p{Hq) cx e ^-'i= 1 0;JJ . 
Using eigenvalue equation UH = AU, p(H, Y\Hq, Yq) can be transformed from matrix space to eigenvalue-eigenvector 
space {A, U} which followed by an integration over ensemble Hq gives 



N 



p(H,Y)cx [I (Afe-A^exp 



k,l;k<l 



N 

-(1/2) J2 A| - (M/2) £ \Xk - A, \ 2 U 2 k U 2 
j=i k<i 



(A2) 



where p = {e^(Y-Y ) _ 

A substitution of eq.(A2) for p in eq.(Al) and using dH = 
{U, A} variables. As eq.(A2) indicates, the behavior of Q™ n . k 



UjUk<i |Afc — Xi\^d\jdUj, gives Q as a function of 
is significantly influenced by the term R = p\Xk — 



Xj\ 2 ElLi \Unk\ 2 \U njl 



present in the exponent of p. Consequently, for a given Y , the dominant contribution to the 
integrals over the variables Uj and Xj in eq.(Al) comes from those regions which lead to R — > 0. Also note that 
the eigenvalue-eigenfunction correlations appear in p only through R. The limit R — > therefore allows a mutually 



independent integration over Xj and U n j variables. As the typical local intensity |C/, 
localization length of the eigenfunction Uk {d as system-dimension), this implies R ~ 

the regions of variable Xj and U n j which contribute to integral depend on mutual competition between p and : 



l k \typical 



|A fc 



A, | 2 



with as the 
Consequently, 



(i) for /x < almost entire region of U n j can contribute to integral (due to < 



|[/ nj | 2 < 1). However only a small 



neighborhood of the order of local mean level spacing, i.e |Afc — Xj\ ^ Dk around A^, contributes to Xj integration. 
Here Dk is the local mean level spacing at eigenvalue Afe. As a consequence, an approximation of repulsion term 
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|Afc — Xj\ ~ D k alongwith the relation X^feLi U n kU* mk = S mn (due to unitary nature of U) gives 

Here r = 0, 1 only, x = 1, -Dfc is the local mean level spacing at eigenvalue X k . 

(ii) for fi > the significant contribution comes from the regions of Xj where \Xj — Ajt | <~ U^C^/^] 1 / 2 . Here 
again, as a typical \U n j\ 2 ~ C~ d < 1) the entire region of E/j can contribute to the integral. Consequently one can 
approximate 

>f (UnjU^Y ( »Y' 2 [5 mn -u n kU* mk y ( 
^^x-W vc fe V (N-iy-iDi [ > 

(One may also consider the contribution from regions where \U n j\ 2 < (p\U nk \ 2 D 2 .) -1 however it is weaker than the 
above). 

By substituting approximations (A3,A4) in eq.(Al), Q™ n . k can be written as (for r = 0, 1 only): 



Qmn-k-X' (N^iy-lps ^Nl{^k,e k ) 



(A5) 



where X = 1 for fj, < and X = M/Cfc for X > /VCk- 

The integral G r (see eq. (30) can also be rewritten in terms of p(H ) and can similarly be approximated: 

G r (x,e)^ ]T J 5% pdH (A6) 

The dominant contribution in this case comes from those regions of integration over Uj and Xj which lead to R = 
p\U n k\ 2 ^2j \Xk — Xj\ 2 \U n j\ 2 present in the exponent of p. (Note, unlike the dominating term R in Q^n-k case : R 



contains only a single component of the k th eigenfunction, namely, U nk , and, the latter takes a fixed value x/\/N.) 
Consequently, for a given Y, G r depends on the mutual competition between \x and x. Reasoning as in the case of 
Qmn;k-> *~*r can be approximated as 

G r « p X v(N -iy- r {N -[xfyP^/D 2 (A7) 
with xo = M _1 for ^\ x \ 2 < 1 an d Xo ~ M 2 for p\x\ 2 > 1. 

APPENDIX B: EFFECT OF MATRIX ELEMENTS PERTURBATIONS ON EIGENVALUES AND 

EIGENFUNCTIONS 

Consider the perturbation of a Hermitian matrix H with matrix elements H k i = ^^^(i)^ 1 Hki-s, eigenvalues 
A„ and cigenfunctions U n , n = 1,2, ..N. By using the eigenvalue equation ^2 m H nm U m j — X n U n j along with the 
ortho- normal condition on eigenvectors i.e. J2j U n jU m j — 5 mn , it can be shown that 



= 2 9ki lu knUin 
On k i- s 

f)TJ . js—1 i 

mp- = — E T — ^-u nm {u* km u ln + (-i)' +1 ur m u ki ). (Bi) 
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The details of the steps used in derivation of eq.(Bl) can be found in [491 ]. 
The set of equations (Bl) can further be used to show following relations: 



and, 



my" Hkl ' s = Xn ( B2 ) 

k,L,s:k<L 

£ - 0, (B3, 



k:l,s:k<l 



(B4) 



E 9kl d 2 U n3 ^nj - , 

2 aff, 2 , ^ (A, - A m ) 2 1 j 

k,l,s;k<l m^j v J my 

■ 9H ^:^ = (B6) 

k,l,s;k<l Kt,S Kl ' S 

E0U n i OUpj U n iU n j . . /T»cr\ 

M , s;fc / fc ' ^^ ;s = ^(A^F (1 ~ ^ (B7) 

E 8U ni dUpj \ - UnmUpm x 

,, -SI *s^SB5:-''Ew^5 1 '" (B8) 
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III. FIGURE CAPTION 

Fig. 1. Distribution P u (u') with u' = [km — (u)]/(ln 2 u) of the local intensity of an eigenfunction near band center 
for AE t (Cubic lattice of linear size L — 13, with hard wall boundary conditions, random hopping and time-reversal 
symmetry) and its BE and PE analog. The parts (a) and (b) of the figure show short and long range behavior of the dis- 
tribution, respectively. The analogs are obtained by the relation I^ 1011 /F a (Q) — I 2 vptcal / Fb(0) — j^ y P tcal J F p (Q) . For 

N = 2197, we find = 0.018, J« 2 = 0.02, and, F a (0) = 0.26, F fc=0 . 02 (0) - (n) 1 ^ , F p=OA (0) = 

0.39. 

Fig. 2 The local intensity distribution for AE nt (Cubic lattice of linear size L = 13, with periodic bound- 
ary conditions, non-random hopping and no time-reversal symmetry) and its BE and PE analogs. In this case, 
jtypical = 0.013,7*^3 =,I^a = 0.00045, and, F o (0) = 0.016, F 6= o.o 3 (0) - W^/ 2 , F p=0 . 4 (0) = for N - 2197. 
Other details are same as in figure 1. 

Fig. 3. Distribution P(I' 2 ) of the rescaled inverse participation ratio I 2 = ln[I 2 /ll yp ] for AE t (same as in Figure 1(a)) 
and its BE and PE analog: (a) short range behavior (lin-lin plot), (b) tail behavior (lin-log plot). Here the BE and 
PE analogs are obtained by the relation Aj ;0 = A/ ; b = A/ ;P . This gives a BE analog of AE different from the figure 1 
although PE analog remains unaffected; the reason lies in almost similar mean level density behavior near band-center 
for the AE, PE cases. 

Fig. 4. Distribution P{I 2 ) of the rescaled inverse participation ratio I 2 for AE„ t (same as in Figure 1(b)) and its BE 
and PE analog. The other details are same as in figure 3. Note the BE analog of AE, PE in this case is different from 
the figure 2. 

Fig. 5. Distribution P{w') of the spatial correlation w' = Inw = \n\zi n z Nn \ 2 between two points belonging to opposite 
end of the sample: (a) short range behavior (lin-lin plot), (b) tail behavior (lin-log plot). The cases compared here 
are AE t (same as in Figure 1(a)) and its BE and PE analog (obtained by the relation A w _ a = A Wt i, = A„,, p ). Again 
the BE analog of AE, PE in this case turns out to be different from the figure 1 but same as in the figure 3. 
Fig. 6. Distribution P{w') of the spatial correlation w' for AE„ t (same as in figure 2) and its BE and PE analog; 
other details are same as in figure 5. Here the BE and PE analogs are obtained by the relation A„, ;0 = A w ^ = A WtP . 
The BE analog of AE, PE in this case is different from the figure 2 but same as in the figure 4. 

Fig. 7. Distribution P(S) of the nearest-neighbor spacing distribution S of the eigenvalues, with (a), (b) showing 
short and long range behavior, respectively, for AE t (same as in Figure 1) and its BE and PE analog. Here the BE 
and PE analogs are obtained by the relation A Eya = A^b = A 6jP . Note the BE analog of AE, PE in this case is different 
from the figure 1 but same as in the figures 3 and 5. 

Fig. 8. Distribution P(S) of the nearest-neighbor spacing distribution S for the case AE„ t (same as in Figure 2) and 
its BE and PE analog (other details same as in figure 7). 
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